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ABSTRACT 


Measurements of mean velocity magnitude and direction as well as 
three-dimensional turbulence intensity were made in the flow over a 
model of an elevated STOL-port. A 1:300 scale model was placed in a 
wind tunnel flow simulating the mean velocity profile and turbulence 
characteristics of atmospheric winds over a typical city environment 
excluding detailed wake structures of possible nearby buildings. Hot- 
wire anemometer measurements of velocity and turbulence were made along 
approach and departure paths of aircraft operating on the runway centerline 
and at specified lateral distances from the centerline. Approach flow 
directions simulated were 0 and 30 degrees to the runway centerline. 
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1. INTRODUCTION 


With the serious consideration of STOL aircraft as a major transportation 
mode, extensive investigations of STOL feasibility have been and are underway. 
One of the important factors in STOL-port design recognized by the FAA [1] is 
the effect of winds on aircraft operating from these facilities. Interest in 
this problem has accelerated in the last few years as evidenced by an 
increased level of investigation. The effects of atmospheric turbulence on 
small aircraft performance has been discussed [2]. Cass et al. [3] 
investigated the characteristics of turbulence about an airport site and 
found agreement on the influence of nearby buildings with a study by Colmer 
[4] of the wake structure behind an isolated hangar. Kurlowski et al. [5] 
have presented an improved turbulence model for use in flight simulations. 

One STOL-port configuration under study is an elevated structure. It 
would have the advantage of getting an aircraft up away rapidly from 
turbulence associated with surrounding structures. One potential problem for 
this type of STOL-port is the influence of the structure itself on the winds 
in the regions through which the aircraft must fly. A preliminary, qualitative 
evaluation of STOL-port design was performed by Parker et al. [6] which showed 
that proper design could substantially alleviate structure induced flow 
problems. The purpose of the present study was to determine quantitatively 
the influence of a proposed elevated STOL-port on the mean flow and turbulence 
over the structure. The results of this study are to be incorporated into a 
piloted-simulator to determine the effect on pilot-vehicle performance. 

The study was performed by placing a model of the proposed elevated 
STOL-port in a wind tunnel specifically designed to simulate atmospheric winds. 
Proper modeling of the atmospheric boundary layer flow approaching the model 
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is essential in providing accurate determination of the flow characteristics . 
Modeling criteria are discussed in Reference 7. Mean velocity and three- 
dimensional turbulence measurements were made along typical approach and 
departure paths for STOL aircraft. Two wind directions were employed: one 

parallel to the flight path and one 30 degrees to the flight path. Measure- 
ments without the model in place allowed the influence of the structure to be 


evaluated. 
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2. EXPERIMENTAL CONFIGURATION 

2. 1 Wind Tunnel 

The wind study was performed in the Environmental Wind Tunnel 
located in the Fluid Dynamics and Diffusion Laboratory at Colorado 
State University, Fig. 1. The tunnel is an open circuit facility driven 
by a variable pitch propeller. The test section is nominally 12 ft 
(3.66 m) wide, 8 ft (2.44 m) high, and 52 ft (15.8 m) long. The roof 
is adjustable to maintain a zero pressure gradient along the test section. 
The mean velocity can be adjusted continuously from 2 to 24 fps (0.61 to 
7.31 mps) with the 10 hp motor and to a maximum velocity of 60 fps 
(18.3 mps) with the 150 hp motor. 

2 . 2 Model 

The elevated STOL-port structure has dimensions of 2200 ft (670.5 m) 
by 1200 ft (365.8 m) by 100 ft (30.5 m) high. In order to provide 
maximum resolution in the data, the largest scale model feasible within 
the physical limitations imposed by the wind tunnel was chosen. Scale 
selected was 1:300. Dimensions of the scaled model were 88 in (2.23 m) 
long, 48 in (1.22 m) wide, and 4 in (10.2 cm) high. The model was placed 
25 ft (7.6 m) from the test section entrance on a turntable to allow ease 
of rotation. The turbulent boundary layer approaching the structure was 
developed with a carefully designed trip at the test section entrance 
12 in (30.5 cm) high consisting of 2.5 in (6.4 cm) diameter tubes 
followed by h in (1.27 cm) gravel uniformly spread over the floor. A 
photograph of the model in-place in the wind tunnel is shown in Fig. 2. 
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No attempt was made to simulate the detailed wake of individual 
structures in the flow approaching the model. Instead, a turbulent 
boundary layer was generated whose characteristics simulated a typical 
atmospheric boundary- layer flow over an urban environment. 
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3. INSTRUMENTATION AND DATA ACQUISITION 

3. 1 Flow Visualization 

Visualization of the flow over the model is helpful in understanding 
and interpreting mean and fluctuating pressures, in defining zones of 
separated and reattached flow, and in identifying size and intensity of 
vortices generated by building corners. Titanium tetrachloride smoke 
was released from sources on and near the model and photographic records 
were made of the flow. Conclusions based on these observations are 
discussed in section 4.1. 

3.2 Mean Velocity and Flow Direction 

To measure three components of mean flow and the mean velocity 
direction a new hot-wire anemometer measurement scheme which has been 
under development at the Fluid Dynamics and Diffusion Laboratory at 
Colorado State University was employed. The procedure utilized a con- 
ventional 45 degree yawed hot-wire sensing element 0.00015 in. (0.0038 mm) 
in diameter with an active length of 0.050 in. (1.25 mm). A Colorado 
State University hot-wire anemometer unit was used with the sensing 
element. The unit is characterized by good stability and high signal 
to noise ratio. A photograph of the probe is shown in Fig. 3. 

Operation of the measurement system required that the hot-wire be 
calibrated for flow speed and angle-of-flow incidence to the wire. A 
Thermo-Systems flow calibrator was used with a Colorado State University 
built probe positioner. The calibration data was fit to a modified 
King's law equation of the form 
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E 2 = A + BU t C f (0) 


f(0) = 1 + a^0 + a^0^ + 


CD 


where and 0 are the velocity and angle to the wire (0 degrees was 

parallel to the wire) of the calibration flow, E is the mean voltage 
output from the anemometer, and A, B, C, a , a^, are calibration 
constants used to fit the equation to the data. Thirty calibration 
points representing 5 velocities at 6 angles were used in a computer 
program to fit the calibration-equation coefficients by least 
squares. Calculations were performed on the Colorado State University 
CDC6400 computer. A typical calibration is shown in Fig. 4. 

To obtain mean velocity data, the probe was mounted in the wind 
tunnel approximately 20 degrees to the horizontal plane and 20 degrees 
to the tunnel vertical centerline plane. Twelve to 15 mean voltages at 
different rotational positions were obtained using an integrating digital 
voltmeter for 30 seconds. This positioned the sensing element at angles 
to the mean flow ranging from about 20° to 70°. Data was reduced by 
finding the three mean velocity components U m , V^, which provided 
the least squares fit of the data to Eq. 1 in a coordinate system 
aligned with the probe axis of rotation. Because the equation cannot 
be solved explicitly for the velocity components, an iteration method 
was devised to converge to the correct solution from initial velocity 
estimates. Convergence was obtained in 6 to 8 interations with simple 
initial estimates. A vector rotation was applied to the velocity compo- 
nents to transform them from the probe oriented coordinate system into 
the coordinate system defined for the STOL-port (Section 3.4 and Figure 5). 
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Horizontal and vertical flow angles were computed from the velocity 
components. All calculations were performed on the Colorado State 
University CDC6400 computer. 

The probe attitude relative to the wind tunnel coordinate system 
was set with an optical alignment system. A laser light source was aligned 
along the data plane being measured. A mirror mounted to the probe 
support reflected the light back to a target placed near the laser. 

Using this system, the probe could be aligned to a predetermined 
attitude within 0.2 degrees. 

The accuracy of the mean velocity measurements was measured in 
such a way that the error included both the variablity inherent in the 
wind-tunnel flow and the uncertainty due to the measurement probe 
itself. The analysis was performed by measuring several vertical pro- 
files in the wind tunnel with the model removed but the roughness 
remaining. The total error was the difference between the assumed zero 
angle of flow and the angle measured with the probe. The standard 
deviation in a (the flow angle in the horizontal plane) was 0.90 
degrees while the standard deviation in 6 (the flow angle in the 
vertical plane) was 0.79 degrees. These figures compare favorably 
to the angular resolution of 5.0 degrees in the horizontal and 1.0 
degrees in the vertical requested in the statement of work defining 
the scope of this investigation. The error in the magnitude of the 
total velocity vector was not measured precisely because the errors 
in the comparative pitot tube measurements appeared from tests to be 
larger than those of the probe itself. The error appeared to be one 


percent or less. 
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3.3 Turbulence Intensities 

Measurement of turbulence intensity was made with a standard 
X-wire probe. The probe was aligned to the mean flow direction deter- 
mined previously with the optical alignment system described above. 

Each of the wires was calibrated individually with the Thermo -Systems 
calibrator. Data reduction used the procedure described by Klatt [12] 
which accounts for differences in wire calibration for each wire. The 

turbulence quantities calculated were u /U , v /U , 

rms vector rms vector 

and w /U where u, v, and w are the three fluctuating velocities 

rms vector b 

in directions X, Y, Z and U vet ^ is the local mean velocity magnitude. 

An independent check of the crossed-wire turbulence results was 
made by measuring at several data locations the longitudinal component 
of velocity fluctuation with a hot-wire placed normal to the flow. The 
results of the comparative measurements are shown by the triangles and 
squares in Fig. 5. The results are in excellent agreement. 

3.4 Data Locations 

Except for vertical profiles to determine the characteristics 
of the approach flow, all data was taken along data lines lying within 
data planes. Two data planes were used: 1) Center Data Plane, a 

vertical plane through the center of the building 600 ft (182.0 m) 
from either side and 2) Side Data Plane, a plane parallel to the 
Center Plane but 150 ft (45.7 m) in from one side. These data planes 
provided information on the conditions an aircraft could experience 
during take-off or landing at the center of the runway and near the 
edge of the structure. 

A data line was established in each data plane along lines 
expected to be used for approach, runway, and departure. The approach 
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segment was defined as a line extending upward at a 7.5 degree angle 
to the horizontal in a downwind direction from a point on the runway 
centerline 350 ft (106.7 m) from the downwind edge of the roof. The 
departure segment was defined as a line extending upwards at a 7.5 
degree angle to the horizontal in an upwind direction from a point on 
the runway centerline 1000 ft (304.8 m) from the downwind edge of the 
roof. The runway segment was defined as a line 15 ft (4.6 m) above 
the runway centerline that connected the approach segment to the 
departure segment. The data lines are shown in Fig. 6. 

Data points were located along the data lines to provide adequate 
resolution of the measurements. Five data points were evenly spaced 
along the runway segment, seven points were located along the approach 
segment and thirteen points, were located along the departure segment for 
a total of twenty-five data points per data line. 

A coordinate system fixed to the STOL-port was defined for data 
reporting purposes (Fig. 6). Y = 0 defined the particular data plane 
of interest. Negative Z represented height above the STOL-port roof 
while positive X was distance into the wind from the downwind edge of 
the roof. This coordinate system is the same as provided in the statement 
of work defining the scope of this investigation. 
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4 . RESULTS 

4 . 1 Flow Visualization 

Structures with sharp corners placed in an atmospheric flow are 
known to generate vortices from leading corners when the wind is at 
an angle to a building side. Increased turbulence in the approach flow 
tends to break down organized vortex structures more rapidly. A primary 
objective of the visualization was to identfy this vortex structure by 
observing its strength and coherence downstream. Smoke flow showed 
significant vorticity with an axis parallel to the structure's long side 
over the model near the windward side for flow approaching at 30 degrees 
to the building long axis. This vorticity emanated from the leading 
corner but was fed along its length by separation of flow along the 
leading edge of the roof. Flow reattached to the roof not far downstream 
from the leading edge, but neither the separation streamline nor 
reattachment region' were well defined. Diffusion of smoke occurred so 
rapidly in the turbulence that the above described phenomena could 
not be adequately captured on film. The influence of this vorticity 
on flight patterns should be confined to the windward side data plane 
for the 30 degree wind. 

Lifting of the flow approaching the structure appeared to extend to 
rather high elevations over the building. A variation in vertical flow 
angle of possible significance to aircraft was observed with the smoke 
for all cases in the region over the structure's leading edge. All 
conclusions obtained from flow visualization were confirmed with the 
velocity measurements reported in sections 4.3 and 4.4. 

4 . 2 Data Scaling 

Velocity data was taken for a reference velocity outside the 
boundary layer of 23.2 fps (7.07 mps) . This provided a wind tunnel 
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Reynolds number ( — ) based on building height and approach velocity 

4 

at that height of 3 x 10 . This is the smallest of the Reynolds 
numbers which can be calculated using building geometry for the length 
scale, and a characteristic approach velocity. At this value of 
Reynolds number, the flow structure should be Reynolds number indepen- 
dent; that is a component of mean velocity at any point in the flow 
normalized by the velocity at the edge of the boundary layer should 
be independent of the approach velocity U . In addition, a turbulence . 
intensity as defined in section 3.3 at any point in the flow should 
be independent of U . Thus the data taken at 23.2 fps (7.07 mps) can 
be used to obtain mean and root-mean- square velocities for any desired 
approach velocity. Mean velocity magnitudes are obtained by simply 
multiplying the nondimensional velocities (i.e. U/U ) by the desired 
free stream velocity. Root-mean-square velocities must, in addition, 
be mutiplied by the nondimensionalized local mean velocity, 

U vector 

— — , since the turbulence quantities were made nondimensional 

oo 

with the local velocity according to standard practice. 

An experimental verification of the Reynolds number independence 
was performed by remeasuring several quantities at a different approach 
velocity. A comparison of mean velocity profiles taken along the upwind 
centerline data line for a 30 degree approach wind for 23.2 fps (7.07 
mps) and 15 fps (4.57 mps) is shown in Fig. 7. The comparison is good. 
The profile shown should be valid for all approach velocities above 
15 fps. A comparison of longitudinal turbulence intensity for a portion 
of the same data line for identical velocities is shown by the circles 
and triangles in Fig. 5. The agreement is within the accuracy required 
for the study. The turbulent intensities are thus valid for the same 
range of approach velocities as the mean velocity. 
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4 . 3 Mean Velocities 

Mean velocity profiles were obtained at the model site with the 

building removed to provide a baseline for comparison of data with the 

model in place. These profiles are shown in Figs. 8 and 9. The 

f Z -Hi 

vertical coordinate shown, - — ^ - , is in the coordinate system 
defined in Section 3.4 and represents the nondimensional distance 
above the ground level. The velocities are shown in both logarithmic 
and semi-logarithmic form to demonstrate the characteristics of the 
boundary layer. For the logarithmic description, the velocity profile 
is expressed in the following form: 

_U_ _ ( Z-H. 0.345 

U 6 ] 

00 

The exponent 0.345 is centered within the range of values expected for 
a city environment. The velocity profile in semi-logarithmic form 
shows a linear profile as is expected for a neutrally stable atmospheric 
boundary layer. The roughness length Z q is 7.2 ft. (2.2 m) in the 
full scale, a very reasonable value for a city roughness. 

Mean velocity magnitudes and directions with the model in place 
are recorded in Tables 1 and 2. Table 1 shows the magnitude of the 
mean velocity vector and its horizontal and vertical angles to the 
reference coordinate system for each data point in the 5 data planes of 
interest. The velocity vector magnitude is recorded as U . /U 
where is the velocity at the edge of the boundary layer. Definitions 

of the horizontal angle a and vertical angle 8 are given in Fig. 10. 
With the model oriented so that the wind vector is from an angle to the 
left for a landing or departing aircraft, an a more positive than the 
prevailing wind angle would indicate a local wind to the right of the 
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prevailing direction. A positive 8 angle indicates an upward 
vertical component in the wind vector. The data planes investigated 
were the center and one side plane for a zero degree approach wind and 
the center and both side planes for a 30 degree wind approaching from 
the left. Table 1 also shows the magnitude of the mean velocity 
vector on the center data plane at zero degree wind azimuth with the 
STOL port removed from the wind tunnel. The a and 8 angles were 
nominally zero for that case. 

A breakdown for the velocity vector into components is shown in . 
Table 2. Data is supplied only for the 30 degree case since the V 
and W components were negligibly small for the 0 degree cases. 

An overall view of the mean flow characteristics along the data 
lines is shown in Fig. 11. The individual arrows represent the 
magnitude and horizontal direction of the mean flow measured at the 
arrow base on three data lines. The strongest effects are seen in the 
upwind side data plane. The effect of the model on the velocity 
magnitude can be seen in Fig. 12. The three curves represent (1) the 
approach flow without the model and (2) with the model, the velocities 
on the approach and departure segments of the center data line for a zero 
degree wind. Reduction in wind velocity with the model in place is 
evident for both cases. 

Of great concern to aircraft is the angle of the flow over the 
structure --in particular, rapid changes in vertical angle of flow may 
cause difficulties for small aircraft. Plots of vertical and horizon- 
tal angles for the mean flow vector are plotted in Figs. 13 through 
16. Figure 13 shows the vertical angle for the data lines in the 
center and side planes for a zero degree wind. Extremum in angle 
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occur at or just off the edge of the structure with little change in 
angle over the center portion of the structure. Figure 14 shows the 
horizontal flow angle for the same conditions as the previous figure. 
No sharp gradients exist in the angles. The approximately one 
degree bias in the center plane data is due to a small cross flow 
tendency in the lower third of the 8 ft (2.44 m) height of the wind 
tunnel flow (the model was aligned precisely with the geometrical 
wind tunnel axis) . 

The vertical flow angle for the center plane and both side planes 
is shown in Fig. 15. All three show the same qualitative features 
with the magnitude of effect differing slightly. Again, the largest 
gradient in angle occurred near the edges of the model. Figure 16 
shows the horizontal angle for the same conditions. The gradient of 
angle is severe only for the windward side plane, an effect also 
noted in Fig. 11. Smoke flow observations tended to show the same 
effect. 

4.4 Turbulence Intensities 

Longitudinal turbulence intensity at the building site without 
the model in place is shown in Fig. 17 as a function of height above 
the surface. The data provides a comparison for data taken with the 
model in place. 

Turbulence intensities in the three directions defined by the 
coordinate system defined in Section 3.4 are listed in Table 3. 
Measurement positions for two data planes at a 0 degree wind and three 
data planes at 30 degree wind were the same as for mean velocities. 
Some characteristics of the turbulent field can be seen from Figs. 18 
and 19. For a 0 degree wind, the longitudinal fluctuations varied 
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smoothly reaching a maximum amplitude of 15 to 17 percent in the two 
data planes at about 80 to 100 ft above the STOL port surface on the 
departure line. The fluctuations varied more for the 30 degree wind 
(Fig. 19) reaching a maximum value of 23 percent near the takeoff point 
on the departure line on the windward data plane. Maximum value on 
the other two planes at 30° were 12 percent. The spatial variability 
of the windward data was confirmed with repeat measurements. The 
turbulence components in the Y and Z directions generally followed the 
same trends as the longitudinal component, but at lower amplitudes. 
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5. CONCLUSIONS . 

Measurements of mean velocity magnitude and direction as well as 
three-dimensional turbulence intensity were made in the flow over a 
modeled STOL port. The flow approaching the model was representative 
of flow over a city environment excluding detailed wake structures of 
possible nearby buildings. The measurements were made along approach 
and departure paths of aircraft operating on the building center line 
and 150 ft in from each side of the structure with wind directions 
parallel and 30 degrees to the flight path. Results of the measurements 
showed the worst environment of those studied to be the flight path 
nearest the windward side of the STOL-port with a 30 degree wind. 

Sizable turbulence intensities and gradients of mean flow angle were 
found, however, even for the centerline with a 0 degree wind. 

The present study did not attempt to modify the flow characteristics 
over the structure with vortex spoiling devices. Relatively simple 
modifications to the local geometry at the building edge which would 
not interfere with aircraft safety should help to alleviate some of 
the adverse flow characteristics seen in this study. Since the elevated 
STOL port was conceived for use near a large city, the wakes of other 
structures in the vicinity may cause larger effects on the flow than the 
structure itself. For this reason, additional wind-tunnel tests with 
realistic city surroundings should be undertaken before suitability of 
the elevated STOL port concept is accepted. 
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MEAN VELOCITY COMPONENTS FOR 30 DEGREE WIND 
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Figure 2. Model installed in the wind tunnel for 30 degree wind. 
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Wire Response, 
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Figure 4. Typical hot-wire calibration 
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Figure 5. Results of hot-wire anemometer turbulence validation tests 






Figure 6. Data lines and coordinate system 
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Figure 7 


Demonstration of Reynolds number independence for velocity 


x (f t ) 




Figure 8. Mean velocity profile in logarithmic form 
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Figure 9. Mean velocity profile in semi- logarithmic form 
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Figure 12. Mean velocity with and without model in place 
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Figure 14. Horizontal flow angles for 0 degree wind 
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Figure 16. Horizontal flow angles for 30 degree wind 
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1 turbulence intensity for 0 degree wind 





Figure 19. Longitudinal turbulence intensity for 30 degree wind 
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